Phase structure of Z2 gauge theories for 
frustrated antiferromagnets in two dimensions 
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In this paper, we study phase structure of Z2 lattice gauge theories that appear as an effective 
field theory describing low-energy properties of frustrated antiferromagnets in two dimensions. Spin 
operators are expressed in terms of Schwinger bosons, and an emergent U(l) gauge symmetry reduces 
to a Z2 gauge symmetry as a result of condensation of a bilinear operator of the Schwinger boson 
describing a short-range spiral order. We investigated the phase structure of the gauge theories 
by means of the Monte-Carlo simulations, and found that there exist three phases, phase with a 
long-range spiral order, a dimer state, and a spin liquid with deconfined spinous. Detailed phase 
structure and properties of phase transitions depend on details of the models. 
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I. INTRODUCTION 

In the last few decades, strongly-correlated electron 
systems are one of the most intensively studied areas in 
the condensed matter physics. One may expect that some 
exotic phase appears as a result of the interplay of strong 
correlations and quantum fluctuations. Concerning to 
the high- Tc cuprates, understanding of the under-doped 
regime is still controversial. Conventional Fermi-liquid 
picture may not hold in that region [1)]. 

Another intensively studied system is quantum mag- 
nets with frustrations. Study of that system has long his- 
tory but its interests recently revived because very inter- 
esting experiments on the new materials like the organic 
Mott insulators k-(ET)2Z (Z=Cu[N(CN)2]C1, etc)[i| and 
X[Pd(dmit)2]2 (X=Me4P, etc)[l, i, ^ have appeared. 
Among them, the insulator with Z=Cu2(CN)3 has no 
long-range order at low temperature 0, Q and it is ex- 
pected that a new type of spin liquid, so called Z2 spin liq- 
uid, is realized there[8|. Another interesting anisotropic 
triangular antiferromagnet is CS2CUCI4. By neutron 
scattering, its spinon-like behaviors were observed [1113. 

To study possibility of exotic states in frustrated anti- 
ferromagnets like the Z2 spin liquid, most studies employ 
the Schwinger-boson representation for quantum spin op- 
erator. As a result, there appear a local U(l) gauge sym- 
metry and also an emergent gauge field. Dynamics of the 
emergent gauge field strongly influences the structure of 
the ground state and low-energy excitations. In the Z2 
spin-liquid scenario, the U(l) gauge symmetry is reduced 
to a Z2 symmetry because of appearance of a short-range 
spin spiral order, and s = | spinous are deconfined and 
appear as a low-energy excitation (llj. In order to obtain 
a conclusive proof of the existence of the Z2 spin-liquid, 
reliable investigation on the gauge dynamics is necessary. 
In the present paper, we shall report results of study on 
the Z2 gauge theories obtained mostly by means of the 
Monte-Carlo (MC) simulations. 

The present paper is organized as follows. In Sec. II, we 
shall introduce models of frustrated antiferromagnets and 



review the Schwinger-boson representation of them. We 
show that their low-energy effective model is a CP^ gauge 
model coupled with an additional doubly-charged vector 
field describing a short-range spiral order. In Sec. Ill, we 
shall show the phase structure of various effective gauge 
models with local Z2 gauge symmetry. To obtain the 
phase diagrams, we calculated "internal energy" , "spe- 
cific heat" , spin correlation functions and instanton den- 
sity by means of MC simulations. There are three phases, 
phase with long-range order, dimer phase, and spin liq- 
uid with deconfined spinous. Section IV is devoted for 
conclusion and discussion. 



II. FRUSTRATED ANTIFERROMAGNTES, 
SCHWINGER BOSON AND EFFECTIVE GAUGE 
THEORY 

A. AF magnets and CP^ gauge field theory 




FIG. 1: Triangular lattice on which the Heisenberg model 
(gHI) is defined. 

Let us start with some specific model of a frustrated 
antiferromagnet on the triangular lattice shown in FiglT] 
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Exchange coupHng in the horizontal bond is J' and the 
others are J. Quantum Hamiltonian Ti is given as 



^L^J^s^■s, + J'Y,s^■SJ 



(2.1) 



and the el- 



where Si in s — ^ spin operator at site i, 
lipsis denotes multi-spin and/or long-range interactions 
between spins, and the other notations are self-evident. 

In the limit J'/ J <C 1, the system reduces to the usual 
antiferromagnets on the square lattice and the ground 
state is expected to have the Neel order, whereas for 
J' I J ^ 1, a new state is expected to appear. In order 
to study the system (|2.ip by field-theory methods, we in- 
troduce the Schwinger boson operators = {a^i,aii) at 
each site z, and then Si is expressed as 



c ^ t 
Si = -a! 



(2.2) 



where a are the Pauli spin matrices. The following local 
constraint must be imposed as the physical-state condi- 
tion in the Schwinger boson Hilbert space. 



(a|,a|, -|-a|^a|,)|Phys) = |Phys) 



(2.3) 



We employ the path-integral methods to investigate 
the quantum system, and introduce CP"'^ variables Zi = 
[z'^i^zyi) — {zii, Z2i) corresponding to a,j, which satisfy 
the constraint 



z^iZ^i + ziiZii = 1, 



(2.4) 



at each site i and is the complex conjugate oi z-]^ia. 
From the Hamiltonian (|2.ip . the partition function is 



given as 



[DzDz](jpi exp 



dr 



(E 



'H(z,z] 



(2.5) 

where r is the imaginary time, Zi = ^ and J[DzDz]qpi 
denotes the integration over CP^ variables z^'s satisfying 
the constraint (|2.4p . TL(z,z) is derived from (|2.ip and 
The above system is obviously invariant under a 

Zi{T) with an 



local gauge transformation Zi^r) 
arbitrary 6'i(r) satisfying 6i{+co) = 0i{—co). 

In the limit J' ^ 0, an effective field theory is obtained 
from the partition function Z in (|2.5p by integrating out 
the high-energy modes of Zi (or ZiS on all odd sites[12,, 
ITsj). The resultant theory is a CP^ gauge model, which 
is described by the following action Sz in the continuum 
spacetime with coordinate a;^t = (a;o = r, xi,a;2). 
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tj,<i' 



where D^^z = (i9^ -I- iA^)z, F^i, = d^Ay - di,A^ with 
emergent gauge field A^. In Eq. ()2.6|) . g, e are couphng 
constants. Bare value of g is independent of the antifer- 
romagnetic (AF) exchange coupling J, but it measures 



the solidity of the AF order, i.e., additional interactions 
that enhance (suppress) the AF order decrease (increase) 
the value of g. On the other hand, the bare value of 
1/e is vanishing for the AF Heisenberg model with only 
the nearest-neighbor (NN) coupling but it acquire a finite 
value due to the renormalization effect of the high-energy 
modes. Multi-spin nonlocal interactions like a ring ex- 
change coupling generate nonvanishing value of l/e[14j. 
Varying the parameters g and e induces a phase transi- 
tion and the structure of the ground state and low-energy 
excitations change drastically through the phase transi- 
tion as we see in the following sections. 

The field theory defined by (|2.6p is obviously invariant 
under a U(l) gauge transformation. The continuum de- 
scription ()2.6p makes it unclear if this U(l) gauge sym- 
metry is compact or noncompact one. As the original 
system of the AF magnets is defined on the lattice and 
transformation parameter 9i{T) is defined mod 2-k, one 
may expect that the model (|2.6p is a compact U(l) gauge 
system, in which topological nontrivial objects like in- 
stantons and vortices can exist. This expectation is qual- 
itatively correct, but contribution from instanton config- 
urations to the partition function is partly suppressed if 
there exists a Berry-phase term, / d?xe^v\d^Fy\ (where 
is the antisymmetric tensor), in the action in addi- 
tion to S'z[ll,[H,[l3l- For the case J' 0, it is not easy to 
calculate the coefficient of the Berry phase, which plays 
a crucial role in the suppression of instantons. We shall 
not consider its effect in the following numerical investi- 
gation, and give comments on it in Sec.IV[18]. 

Phase structure of the CP''^"^ field theory has been 
studied by the 1/iV-expansion and numerical methods [l9l. 
[20l . [2]| . For the compact U(l) gauge case, a lattice- 
regularized version of (|2.6p is quite useful for investiga- 
tion on the CP^ gauge model, and its action is given as 
follows, 



Cl 

2 



C2 



H-c.c, 



^ ^ Ux , fi Ux -\-fi,}y Ux -{-v.^ Ux , 1/ 

(2.7) 



where x denotes site of the cubic lattice, and the cou- 
pling Cl corresponds to l/g"^ and C2 to Phase 
diagram has been obtained in the ci — C2 plane. See 
FigEl There are two phases separated by the critical 
line Cl ~ Cic(c2), one of which corresponds to the Neel 
state for ci > Cic(c2) and the other state is a dimer state 
Cl < Cic(c2) in which the spinon Zx is confined to a spin- 
triplet excitation z^j-az^. The phase transition across the 
transition line is of second order, and it belongs to the 
universality class of the 0(3) nonlinear sigma model in 
three dimensions (3D) (for small to medium values of 
(2.6) ^'^)- The s = ^ AF Heisenberg model corresponds to 
Cl > cic, and the ground state has the AF long-range 
order. By introducing an inhomogeneity in the exchange 
coupling J that enhances dimerization, the value of ci 
in the effective model (|2.7p is decreased and the phase 
transition takes place from the Neel to dimer states(23|. 
Recently, numerical study on the inhomogeneous SU(2) 
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AF Heisenberg model, which is essentially the same with 
that studied in Ref. [22!] . was performed quite in detail 
and the existence of the phase transition from the Neel 
to dimer states was verified |23l|. Phase transition belongs 
to the universality class of the 31? 0(3) nonlinear-sigma 
model, as predicted by the study of the effective lattice 
model 
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of the effective field theory and we expect that details of 
the lattice model does not influence substantially physi- 
cal results like phase structure and critical behaviors by 
the unversality-class argument. For the quantum SU(2) 
AF magnets, it is known that the results obtained by 
the effective CP^ lattice model (|2.7p are in good agree- 
ment with those obtained for the original AF Heisenberg 
model, as we explained above. Furthermore, phase struc- 
ture of the lattice CP" (n — 1 , • • • ,4) models obtained 
by the MC simulations is the same with that obtained 
by the 1/A^-expansion for the CP^~^ field theory in the 
continuum spacetime[l^ These facts encourage us 
to apply the same methods to more complicated quantum 
spin systems like triangular AF spin systems with frus- 
trations. More comments on the reliablity of the methods 
will be given in Sec. IV, after showing the main results of 
the present study in the following sections. 

B. Effect of frustrations 



FIG. 2: Phase diagram of the U(l) gauge theory of CP^ 
spinons[20l|. There are two phases. 

As shown in Figl^l the deconfined Coulomb phase does 
not exist in the model (|2.7p of the U(l) gauge theory. Ap- 
pearance of the Coulomb phase requires long-range and 
nonlocal interaction of gauge field C/x^, which may be 
generated by the coupling with gapless fermions[2^. l25t . 
In the pure quantum spin models without doping of holes, 
the deconfined phase is expected to appear by introduc- 
ing frustrations because the Higgs mechanism is expected 
to take place by the appearance of the (short-range) spi- 
ral order. In that case, the U{1) gauge symmetry spon- 
taneously breaks down to Z2. It is known that the de- 
confined phase exists in the 3D Z2 gauge models. There 
are interesting studies on spin liquids with deconfined 
spinous in the framework of the Z2 gauge model. How- 
ever, detailed and reliable study on the phase structure 
of the Z2 gauge models relevant to the frustrated spin 
systems is still lacking. We study this problem in this 
paper. 

Before going into details of the study on the frus- 
trated AF magnets, let us comment on the validity of 
the present methods using the lattice field theory for 
studying AF magnets. To define quantum many-body 
systems without ambiguities, an ultra-violet (UV) reg- 
ularization is necessary. In quantum spin models like 
(|2.ip . the spatial lattice naturally gives such an UV reg- 
ularization. In the present approach, we first study the 
original model carefully and identify the relevant modes 
in the low-energy and low-momentum region. Through 
these observations, we obtain an effective field theory in 
the continuum spacetime. Then in order to study the 
effective field theory nonperturbatively (e.g. by means 
of the MC simulations), we reformulate it by using a 
spacetime lattice as a systematic regularization. Struc- 
ture of the lattice model is deterimined by the symmetry 



The effect of the frustration in the AF magnets (|2.ip 
can be studied in the framework of the CP^ gauge field 
theor y w hose action has the following term Sa in addition 
to 5,|lH, 



Sa = 



J 1 n ^9a 



a=l,2 
-iAaZdaZ + C.C. 



toaIAcI 



A|A„|4 
(2.8) 



where A^ {a =^ 1,2) is a doubly- charged spatial vector 



field, D^if^ 



+ 2iA^, and Za{x) = eabZb{x) (ei 



—621 — 1, £11 = £22 = 0). Origin of the new term Sa 
is as follows. The J'-term in Ea. (|2.ip generates terms 



like ^J2ij ■ %l the effective field theory, where 
the extra factor 1/J comes from the redefinition of the 
imaginary time t — > t x (a J) (a =lattice spacing=often 
set unity). After inserting the following identity into the 
path-integral representation of the partition function, 

rfAydAy e-4(:FA..-'^.-^=.)(:f A..+»2.-^.) = constant 

(2.9) 

the above quartic term of Zi is decoupled by a Hubbard- 
Stratonovich field Aq,. By the effects of renormalization 
of high- momentum modes, the extra terms in Sa and 
renormalization of the mass, which preserve the local 
U(l) gauge symmetry, appear for describing low-energy 
behavior of the svstem|26|. 

Physical meaning of Sa becomes transparent by con- 
sidering the case toa < 0. In this case, we expect the non- 
vanishing expectation value of the field A^,, i.e., (Aq) ^ 0. 
By solving the field equation derived from the action 
Sz+Sa, it is straightforward to verify that the low-energy 
configurations arc given by. 



:'^^'<^>-" + eabt^,(x)e-*f'<^>-"~V (2.10) 
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where Va{x) (a = 1,2) is a slowly varying complex field 
satisfying [^^(a;)^ = l[23,[2l|. For the configurations 

given by (PTO)) . the SU(2) spin field S{x) = z{x)az{x) 
has the following form, 

S{x) = fii cos(2g^(A) ■x)+n2 sin(2g2(A) • x), 

ni = Re[{)(7w], n2 — Im[i;(Tu], 

nl = nl = 1, ni • = 0. (2.11) 

On the other hand, the "spin-nematic field" 713 = ni x n2 
is given as 713 = vav. It is obvious that S{x) in (|2.1ip 
corresponds to a spiral state if (ni) 7^ 0, {712) 7^ 0. 

By substituting Eq. (|2.10p and Aqq, — (Aq,) into the 
continuum action 5*^ + ^a, low-energy effective theory 
is obtained. Condensation of Aq. not only generates the 
spiral state of S but also a finite mass of the gauge field 
A^. CP^^^ model with a massive "gauge field" has 
been studied in the continuum spacetime by the 1/A^- 
expansion, but the obtained results are not reliable for 
the case of finite N (in particular the case N = 2) because 
an important effect at 0{1/N) coming from topological 
excitations is totally ignored there In fact, the con- 
densation of Aq preserves the local Z2 gauge invariance 
of the system because it carries double charge, and there- 
fore the topological nontrivial excitation carrying a half- 
magnetic quantum, dubbed vison, exists as a low-energy 
excitation 30]. Also in Ref. 31], a quantum phase transi- 
tion between a spin liquid with deconfined spinous and 
magnetically orderd state was studied, and various phys- 
ical quantities were calculated by the l/7V-expansion in 
an effective CP^~^ field theory with a global U(l) sym- 
metry. There it is assumed that effect of the vison can be 
ignored. Our study of the gauge model with the full Z2 
gauge symmetry in the present paper will show that this 
assumption is correct. See, for example, the calculation 
of the instanton density in Sec. III. 

In the rest of the present paper, we shall study the 
effective Z2 gauge theories obtained by substituting A = 
(A) and Eq. fTTU)) into the action -h ^a- To this end, 
we reformulate it by using the lattice regularization that 
preserves the local Z2 gauge symmetry. We use a cu- 
bic spacetime lattice because frustrations coming from 
AF coupling on the triangular lattice has disappeared by 
using the parameterization (|2.10[) . The resultant lattice 
model is explicitly given by the following action, 



(I2.12P is given as 



^(C3) 



Cl 

2 



C-2 \ ^ 



X,fJ,<U 



(2.12) 



where we explicitly show the dependence of the parame- 
ter C3 in A{c3), as we study the model with fixed values 
of C3 in the following section. From the above consider- 
ation, C3 cx (A)^. Partition function of the gauge model 



^Gaugc= J[DzDz]cpi[DUDU] expA(c3). (2.13) 



It is obvious that the system (|2.12p has a local Z2 gauge 
symmetry instead of the U(l) symmetry. Then we call Vx 
ZaCP^ boson. In the limit C3 — > 00, configurations of the 
gauge field are restricted to Ux,fj. = ±1 and the model re- 
duces to a Z2 gauge system. In the case C3 — > 00, ci = 0, 
the system is the pure Z2 gauge model in 3D, which is 
dual to the 3D Ising model and has a second-order phase 
transition from the confined to deconfined "Coulomb" 
phases as C2 is increased. This is in sharp contrast to the 
U(l) gauge model in 31?, in which only the confined phase 
exists. As the deconfined phase corresponds to spin liq- 
uid with weakly interacting spinous, one may expect re- 
alization of a fractionalization phenomenon in frustrated 
AF magnets. In the following sections, we shall study 
phase structure of the model (|2.12p by means of the MC 
simulations. 



III. NUMERICAL STUDIES 
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FIG. 3: Phase diagram of tlie Z2 gauge theory of CP^ spinons. 
There are three phases. 



A. Z2 lattice gauge model of Z2CP^ spinon 

We first study the Z2 gauge model coupled with the 
field Vx that corresponds to the limit C3 — > 00 of the 
model (|2.12p . It is known that Z2 gauge model coupled 
with single component Higgs boson describes nematic 
phase transition, and its phase structure was studied 
by both analytical and numerical methods [s^]. In these 
studies, importance of topological line defects (world 
lines of vison) was emphasized. 

In order to investigate the phase structure of the Z2 
gauge model, we defined the model on the cubic lattice of 
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FIG. 4: E for C2 = 0.5. There is sharp discontinuity at ci ~ 
1.35 that indicates a first-order phase transition. System size 
is L = 24. 



size with the periodic boundary condition and calcu- 
lated the "internal energy" E = {A{oo))/L^, the "specific 
heat" C = ((A(oo) - Ef)/L^, etc. We used the stan- 
dard Metropolis algorithm for the MC simulations [s^. 
The typical statistics used was 10^ MC steps for each 
sample, and the averages and errors were estimated over 
10 ~ 20 samples. Average acceptance probability was 
about 40 ~ 50%. 

The obtained phase diagram is shown in FiglSl There 
are three phases, and calculation of various physical 
quantities gives the following identifications. 

1. In the phase I, there is no AF long-range or- 
der and the gauge dynamics is realized in the 
confined phase. Low-energy excitations are spin- 
triplet bound states of the spinor Vx (triplon), i.e., 
ni{x), n2{x) in Ea. (|2.1ip . We call this phase tilted 
dimer state. 

2. In the phase II, there exists the magnetic long- 
range order of Vx, which corresponds to the spiral 
order of S{x), i.e., (ni) ^ 0, (^2) ^ 0. The gauge 
dynamics is in the Higgs phase because of the con- 
densation of Vx ■ Low-energy excitations are gapless 
spin wave described by uncondensed component of 

Vx- 

3. Phase III represents the paramagnetic spin liq- 
uid state. As for gauge dynamics, a deconfined 
"Coulomb phase" is realized, and the number of 
topological vortices is conserved. Low-energy exci- 
tations are massive spinon Vx- 

We first show the numerical calculations of E and C for 
establishing the above phase diagram in FiglJl We first 
focus on the transition from phase I to II. In FigHl we 
show calculation of as a function of ci with C2 = 0.5. 
It is obvious that there exists a sharp discontinuity at 
Cl ~ 1.35, which indicates a first-order phase transition. 
In order to verify it, we measured distribution of values 
of A{oo) generated in the MC steps, Af[oo;E], which is 
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FIG. 5: Distribution of A{oo) for C2 = 0.5 and ci close to 
the phase transition point. The double-peak structure at 
Cl = 1.3545 confirms the existence of the first-order phase 
transition. System size L = 24. 



defined as 

^Gaugc = [ dE [ [DzDz]cpi[DUDU] e^(°°)(5(A(oo) - i;) 

(3.14) 



= J dEJV[oo;E]. 



We show the result near the critical point in FiglS] 
At Cl = 1.3545, JV[oo;E] has a double-peak structure, 
whereas the others have a single peak. From this result, 
we judge that the first-order phase transition takes place 
at Cl = 1.3545. 

We also measured the instanton density p{x) for C2 — 
0.5 as a function of ci. p{x) is defined as follows for 
the gauge field configuration Ux,f_i — e^^^-'^, 9x.fi ~ 
0, 7r[20|, [341 • First we consider the magnetic flux Ox.^iy 
penetrating plaquette {x,x + p,,x + p + v,x -\- v) 



^X.IIU 



(-27r < e,,^, < 27r). (3.15) 



We decompose Qx.p.f into its integer part Ux.^u, which 
represents the Dirac string (vortex line), and the remain- 
ing part Qx.fiu, 

'dx.fj.u = '2'!mx,fi„ + Qx.fj^^, (-TT < 6x,p;y < tt). (3.16) 
Then instanton density p{x) at the cube around the site 
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FIG. 6: Instanton density for C2 = 0.5 as a function of ci. At 
the phase transition point ci ~ 1.35, it changes its behavior. 
System size L = 24. 

x+2-|-i + |of the dual lattice is defined as 

1 . , (3-17) 

where is the antisymmetric tensor. From the above 
definition, it is obvious that p = (|p(2;)|) measures proba- 
bility of creation/annihilation of magnetic vortex. In iD 
Z2 gauge theory, magnetic vortices in iD can be regarded 
as world Hues of flux quanta dubbed vison. Nonvanish- 
ing value of p means that the number of visons is not 
conserved, and therefore condensation of the vison. The 
result of calculation of p is shown in Fig[6l There is a 
sharp discontinuity at the phase transition ci ~ 1.35. In 
phase I, finite value of p means large fluctuations of the 
gauge field and spinon Vx is confined to gauge-invariant 
composites, {i = 1,2,3). This phenomenon is some- 
times called dual Meissner effect. On the other hand in 
phase II, p is strongly suppressed and the topological or- 
der exists. Later study on the spin correlation function 
reveals that this suppression is due to Higgs mechanism 
by the condensation of Vx. 



-1.8 
-1.9 

-2 
-2.1 
-2.2 
-2.3 
-2.4 











— size24 






H3- size32 






— size40 

















0.7 0.72 0.74 0.76 0.78 0.8 
cl 



FIG. 7: E for ci = 0.3 as a function of C2. System size 
L = 24, 32, 40. There is almost no system-size dependence. 



FIG. 8: C for ci = 0.3 as a function of C2 with L = 24, 32, 40. 
Its system-size dependence indicates that the phase transition 
is of second order. In addition to the standard MC simu- 
lations, we used multi-histogram methods to obtain reliable 
values of C near the phase transition point [s^. 
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FIG. 9: FSS for ci = 0.3. AU data for L = 24, 32, 40 can be 
fit by single function 0(a:). 

Next we consider the phase transition from phase I to 
III. We show E and C for ci = 0.3 in Figsl7]and[Hl The 
results indicate that there exists a second-order phase 
transition at C2 — 0.76. By the finite-size scaling (FSS) 
hypothesis for C, 

CL{e) = L'''''cp{L^'''e), (3.18) 

where Cl is the "specific heat" of system size L, and 
e = (c2 — C2oo)/c2oo with C200 (the critical coupling for 
L — !■ 00), we estimated the critical exponents v, a by us- 
ing the FSS ((XTg)) and obtained v = 0.63, cr = 0.17 and 
the critical coupling C200 = 0.76. The obtained scaling 
function 4){x) is shown in FiglHl These values are very 
close to those of the pure Z2 gauge model that are ob- 
tained from the data of the iD Ising model by duality [s^. 

We also measured instanton density p and show the 
result in Fig llOl p is a decreasing function of C2 and 
changes its behavior at the phase transition point C2 ~ 
0.76. 

Finally, let us consider the phase transition from the 
phases II to III. Obtained E has no system-size depen- 
dence. System-size dependence of C is shown in Fig llll 
from which we judge that the phase transition is of sec- 
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FIG. 10: Instanton density for ci — 0.3 as a function of C2. At 
the phase transition point C2 — 0.76, it changes its behavior. 
System size L = 24. 



FIG. 12: Correlation function of the spin field ni{x). It has 
a long-range order only in phase II. 



ond order. By the FSS, the critical exponents are esti- 
mated as v — 0.65, a = 0.156, C2oo = 0.93. This value 
of v should be compared with that of the 0(4) nonlinear 
sigma model in 3D, 1^0(4) = 0.75. At present, it is not 
clear for us if the above two phase transitions belong to 
the same universality class. 
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FIG. 13: Correlation function of the spin field 722(3;). It has 
a long-range order only in phase II. 



studies[20|, we employ a Fourier transformation of the 
spin field, e.g., ^3(2;), 



FIG. 11; C for C2 — 1.5 as a function of ci with L = 
16, 24, 32, 40. Its system-size dependence indicates that 
the phase transition is of second order. 

In order to see (non) existence of the magnetic long- 
range order (LRO), we measured correlation functions of 
the spins ni{x), n2{x) and n3{x). They are defined as, 

G.,{r)^^J2{n,{x + r)-n,{x)), i = 1,2,3. (3.19) 

X 

We exhibit the results in Figs fT2l [T3l and [T4l It is obvious 
that only in phase II, they have the LRO. This LRO 
indicates a nonvanishing expectation value of Vx, {vx) ^ 
0, in phase II. This understanding is supported by the 
measure of p, which shows that the gauge dynamics is in 
the Higgs phase in phase II. 

We also calculated the spin gap as a function ci for 
C2 = 1.5, i.e., from the spin liquid to spiral state. It 
is difhcult to estimate the spin gap directly from the 
spin correlation functions Gi(r). Then as in the previous 



Mxo;pi,P2) = e'P'^'^'P'^'M^)- (3.20) 

Xl,X2 

In the continuum limit, the correlator of n3(xo;pi,p2) 
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FIG. 14: Correlation function of the spin-nematic field 713(2;). 
It has a long-range order only in phase II. 
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FIG. 15: Spin gap as a function of ci for C2 = 1.5. 



FIG. 16: Phase diagram of the U(l) gauge theory of CP^ 
spinous with C3 = 0. There are two phases. 



behaves as 



{n3{xo;pi,p2) ■ m{Q;pi,p2)) 

where p^ = X]i=i2 3^'i- practical calculation 

on the lattice, we put pi — P2 = 27r/L, and mea- 
sured a/pi + P2 + from the correlation function of 
n3{xo;pi,P2)- We show the result Mg as a function of 
ci in Fig[T21 It is obvious that the spin gap Mg is a 
continuous decreasing function of ci and is vanishing for 
ci > cic ~ 1. This result means that the spin excitation 
has a finite gap in the Z2 spin liquid, whereas the spin 
wave in the spiral state is gapless. 

From all the above calculations, we obtain the phase 
diagram shown in FigO 



B. U(l) gauge field coupled to Z2CP^ spinon: 
cs — case 

Let us consider the case C3 = of the system (|2.12l) . 
The system with ci = is nothing but the pure compact 
U(l) gauge model in 31?. It is well-known that there 
is no phase transition and the system is always in the 
confined phase, though there is a crossover from dense- 
instanton to dilute-instanton regimes as the parameter 
C2 is increased. 

We show the obtained phase diagram of the system 
(|2.12p with C3 = in FiglHl There are two phases, (i) 
phase I is the dimer phase with confinement of spinous 
and without any long-range order, (ii) phase II is the spi- 
ral state with the condensation of Vx- Phase transitions 
separating these two phases are of first order as the cal- 
culations of E in FigsfTTl and [TBI indicates. In order to 
verify this observation, we measured the distribution of 
A{0), Af[0]E], in the MC steps. From FigdU it is ob- 
vious that on the critical line Af[0; E] has a double-peak 
structure whereas it does not off the critical line. 

The second-order phase transition that exists in the Z2 
gauge model of spinous from the confined to deconfined 
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FIG. 17: E for C2 — 0.4. There is sharp discontinuity at 
Cl — 1.44 that indicates a first-order phase transition. System 
size is L = 24. 



phases disappears in the system with C3 = 0. There is a 
crossover line emanating from the crossover point of the 
pure compact U(l) gauge model in 3D. See Figl^Ol This 
fact strongly influences structure of the ground state and 
low-energy excitations of the original spin system. Ab- 
sence of the deconfined phase means that the spin liquid 
phase does not exist in the present case. We measured 
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FIG. 18: E for C2 — 1.8. There is sharp discontinuity at 
Cl — 0.95 that indicates a first-order phase transition. System 
size is L = 24. 
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FIG. 19: Distribution of A{0) for C2 = 1.8 and ci close to the 
phase transition point. The double-peak structure at ci = 
0.947 confirms the existence of the first-order phase transition. 
System size L = 24. 
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FIG. 20: C for ci = 0.5. System size L = 12, 18, 24. There is 
no system-size dependence, i.e., the second-order phase tran- 
sition in the Z2 gauge model reduces to a crossover. The de- 
confined spin-liquid phase does not exists in the model with 
C3 = 0. 



the instanton densityp to verify the above conclusion. In 
particular in the dilute-instanton regime C2 > 1.5, the 
instanton density p is small even in the confined phase 
ci < Cic(c2), but it decreases rapidly at the phase tran- 
sition ci = Cic(c2) to the Higgs phase. See FigEH On 
the other hand for ci — 0.5, p is a decreasing function of 
C2 but does not exhibit any anomalous behavior at the 
crossover C2 — 1.5. See Fig|22l 

We calculated the spin correlation functions in each 
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FIG. 21: Instanton density for C2 = 1.8 as a function of ci. 
At the phase transition point ci ~ 0.95, there is a sharp 
discontinuity. System size L — 24. 
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FIG. 22: Instanton density for ci = 0.5 as a function of C2. 
There is no anomalous behavior at the crossover C2 — 1.5. 
System size L = 24. 



phase and verified that the spin LRO exists only in the 
Higgs phase ci > cic- 



C. Massive U(l) gauge field coupled witli Z2CP^ 
spinon: C3 = 1.0 case 
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FIG. 23: Phase diagram of the U(l) gauge theory of CP^ 
spinous with C3 = 1.0. There are three phases as in the Z2 
gauge model. 
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FIG. 24: C for ci = 0.3. System size L = 12, 18, 24. 



Finally let us consider the case C3 = 1.0. We also in- 
vestigated the case C3 — 2.0 and obtained similar results. 
In the present case C3 > 0, the gauge field Ux,fj. is a U{1) 
variable, but local U(l) gauge symmetry is explicitly bro- 
ken down to Z2 by both the hopping term of and the 
mass term of the gauge field, i.e., U^^^+c.c. It is expected 
that the mass term of the gauge field is a relevant per- 
turbation and therefore the phase structure of the system 
C3 > is qualitatively the same with that of the Z2 gauge 
theory studied in Sec. III. A. 

We show the obtained phase diagram in Figl23l There 
are three phases similarly to the Z2 gauge theory of 
spinous, as it is expected. In FigslSH and [25l we show 
C as a function of C2 for ci = 0.3 and the FSS scaHng 
function obtained from these data. Critical exponents 
are estimated a.s 1/ — 1.26, a = 0.43 and C200 = 1-44. 
From the above result, we think that the present phase 
transition does not belong to the universality class of the 
3D Ising model. At phase transitions from the tilted- 
dimer to spiral phases, E, the distribution of Af[l; E] and 
the instanton density p have similar behaviors to those 
in the previous cases of the first-order phase transition. 
In Figsl26] and [27l we show the result of the instanton 
density. 
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FIG. 25: FSS for ci = 0.3. All data for L = 12, 18, 24 can 
be fit by single function </)(a;). 



FIG. 26: Instanton density for C2 = 0.5 as a function of ci. 
At Cl ~ 1.1, there is sharp discontinuity corresponding to the 
first-order phase transition. System size L = 24. 
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FIG. 27: Instanton density for ci = 0.3 as a function of C2. 
At C2 — 1.4, p changes its behavior. System size L — 24. 



IV. CONCLUSION AND DISCUSSION 

In the present paper, we derived effective gauge models 
that describe low-energy properties of antiferromagnets 
with frustrations in 2D, and studied their phase struc- 
ture mostly by means of MC simulations. We found that 
generally there are three phases in the models, (i) phase 
of the tilted dimer state with spin-triplet excitations, (ii) 
the spiral state with gapless spin wave, (iii) the spin liquid 
with weakly interacting spinous. We identified the order 
of the phase transitions and estimated values of the criti- 
cal exponents of the second-order phase transitions. The 
investigation suggests that for the spin liquid to appear, 
multi-spin and nonlocal interactions are necessary in the 
original spin systems. 

In order to verify the validity of the above results, it 
is important and also possible to study spin systems on 
layered 3-D triangular lattice at finite temperature (T) 
by means of the Schwinger boson (CP^) representations. 
In this case, the systems can be studied directly with the 
spatial lattice as a regularization. In the path-integral 
representation of the partition function Z in (|2.5[) at fi- 
nite T, the r-dependence of Zi is ignored. Then the 
path-integral over CP^ variables Zi's in Z can be per- 
formed without any difficulties by the MC simulations. 
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At present we are studying these systems, and have ob- 
tained prehminary results that support the conclusion in 
the present paper [st!]. 

In the present paper, we mostly focus on the (short- 
range) spiral state with (A^) cx {zi ■ Zj) = Aq ^ 0. 
There is another possibility of canted state like (A^ j^j) = 
(— )*Ao ^ 0. This state can be regarded as a state with 
a ferromagnetic order in the AF background. This state 
also breaks the U(l) gauge invariance down to the Z2 
as the (short-range) spiral state does, and therefore re- 
sults obtained in the present paper are expected to be 
applicable to the canted state. 

Finally let us comment on effects of the Berry phase. 
As we explained in Sec. II, the Berry phase appears after 
integrating out the high-energy modes in the path inte- 
gral in order to derive the effective field theory of AF 
magnets. The Berry phase may play an important role 
though qualitative phase structure is not changed by its 
existence ^38j. Whether the suppression of the instantons 
occurs by the Berry phase strongly depends on its co- 
efficient. For example in the inhomogeneous AF Heisen- 
berg model on a square lattice, the coefficient depends on 
the magnitude of the inhomogeneity and is generally an 



irrationalf^^l. Suppression of instantons does not occur 
in that case and the Neel-dimer phase transition belongs 
to the universality class of the classical 3D 0(3) nonlin- 
ear sigma model, which is equivalent to the CP^ gauge 
model (12. 6p without the Berry phase. This result was 
verified by the numerical study of the inhomogeneous 
SU(2) AF Hciscnbcrg model. We expect that nonvanish- 
ing frustration coupling J' gives a similar effect on the 
Berry phase's coefficient because the most dominant NN 
spin pair configuration is shifted from Zia = tabZjb ihj = 
site, a,b ~ spinor indices) on path-integrating out high- 
energy modes. If this is the case, the Berry phase gives 
only negligible effects on critical behavior of the systems 
under study, and the FSS used in the present paper gives 
reliable estimation of the critical exponents [4ll|. 
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